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(57) Abstract 

The present invention of interfering continuous-wave signals are automically detected in the demodulated output of a quadrature 
amplitude modulation data communications receiver. The transmitter (16) is re-tuned in order to avoid the detected interferer. Statistical 
information concerning the received signal is accumulated at a receiver (18) and is analyzed to determine whether a reduction in the 
received bit-error-rate is due to the presence of an interferer. If an inteferer is determined to be present within the signal transmission 
bandwidth, then the potentially available transmission band is spectrum analyzed to find a new portion of that band in which the signal can 
be transmitted free of the interferer. The results of the statistical analysis and the spectral analysis are conveyed back to the transmitter 
(16) which then readjusts the transmission parameters and forwards the new parameters to the receiver (18). The system resumes operation 
in the new transmission band with no interference and with the bit-error-rate restored to its required level. 
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INTERFERENCE DETECTION AND AVOIDANCE TECHNIQUE 

Technical Field 

The present invention relates generally to data 
5 communications, and more particularly, to a method and 
apparatus for automatically detecting the presence of an 
interfering signal in a data communications receiver and 
retuning a transmitter to avoid the detected interfering 
signal. 

10 

Background Art 

To meet the growing need for ever increasing digital 
data bandwidth for new subscriber services (e.g., high- 
data-rate internet service, video telephony, and high 

15 definition television) telephone companies are looking 
towards the use of very-high-speed digital subscriber 
lines (VDSL) . VDSL provides a means to carry such data 
into the home or business over the existing copper wires 
used by plain old telephone service (POTS) . VDSL 

20 carries digital data at bit rates reaching 52 Mbps using 
carrierless AM/ PM (CAP) modulation, which is related to 
quadrature amplitude modulation (QAM) . Modulation 
contellations can reach 256 symbol points with symbol 
rates up to 648 0 Kbaud. The CAP modulation 

25 constellations look the same as the corresponding 

baseband QAM constellations. The signals are carried 
within the spectral band up to 30 MHz, which is above 
the frequency band used by POTS on the same wire pair. 
Due to the high loss of the copper telephone wires at 

30 these frequencies, the VDSL signals are carried on the 
telephone wires only over the "last mile" (i.e., the 
last segment of copper wire between a central office and 
the user's premises) . 
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Within the spectral band of 0.5 MHz to 3 0 MHz, the 
VDSL transmissions over the telephone wires are subject 
to radio frequency interference (RFI) from AM radio 
broadcasts and amateur radio transmissions, which can be 
5 picked up by the telephone wires acting as receiving 

antenna. This interference can corrupt the received and 
demodulated CAP/QAM signals, causing symbol and bit 
errors in the received data. The AM broadcast signals 
remain indefinitely at the same frequencies in any 
10 neighborhood within range of the local AM radio 

transmitters. The amateur radio signals, however, will 
vary in frequency throughout the day as the operator 
tunes his transmitter to optimize the HF-band skywave 
propagation conditions over the communication links he 
15 is attempting to maintain. There are several narrow 

spectral bands allocated to amateur radio within the 0.5 
MHz to 3 0 MHz band carried by VDSL transmission lines. 
A VDSL transmission line will normally be interfered 
with only by amateur radio transmitters within a few 
20 hundred feet of the transmission line. 

Narrowband interferers, such as AM radio broadcasts 
and amateur radio signals, are not the only interferers 
that corrupt the CAP/QAM signals. Wideband noise that 
extends across a large portion of the VDSL spectrum may 
25 also interfere, causing undesirable changes in the 
carrier-to-noise ratio (CNR) and the bit error rate 
(BER) . Wideband noise can be caused, for example, by 
electrical machinery, internal combustion engines (e.g., 
lawn mowers) or fluorescent lights operating in the 
30 immediate vicinity of a cable drop to the premises. 

Narrowband interferers tend to be more predictable and 
hence somewhat easier to avoid, for example, by 
detecting the frequency of the interferer and subsequent 
adjustment of the CAP/QAM VDSL signal for avoidance 
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thereof. Wideband interferers, on the other hand, are 
somewhat more random, and cannot generally be accounted 
for prior to data transmission and are therefore more 
difficult to avoid. Especially difficult is detection 

5 and characterization of the interferers without 
interruption of the data signal for analysis. 

Thus, what is needed is a method and apparatus for 
recognizing the difference between narrowband 
interferers and wideband noise sources without relying 

10 upon spectrum analysis of the received signal except to 
locate a new transmission band. This is accomplished by 
the present invention whereby a determination is made of 
the presence or absence of a narrowband or continuous 
wave (CW) interferer, and in many cases inferring its 

15 carrier frequency without the need for interrupting the 
transmitted data stream. The transmitted data stream is 
only interrupted when an inteferer is detected and a 
spectrum analysis is needed to find a new transmission 
band. The interferer presence determination can 

20 generally be accomplished when the interferer level is 

low enough that the presence of the interferer cannot be 
inferred by visual examination of the constellation plot 
of the demodulated signal. When a narrowband signal is 
present, the shape of the symbol clusters in the 

25 constellation plot changes from a centrally peaked 

distribution, characteristic of a signal in Gaussian 
noise, to a ring as shown in FIGs . 1A and IB. The 
present invention recognizes this transition 
automatically and at low interference levels where the 

30 transistion to a ring is not yet visually apparent. 

Adjustment of the transmitted symbol rate, bandwidth and 
carrier frequency is provided in order to avoid the 
troublesome interferer and other interferers that may be 
present in the signal environment. 
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Disclosure of Invention 

According to a first embodiment of the present 
invention, a method of non-intrusively detecting whether 
5 an interference signal is present with a received signal 
of interest is provided. The interference signal is 
detected by accumulating statistical information 
associated with the signal of interest. The statistical 
information is analyzed to determine whether an 
10 interference signal is present and is intrusive. If an 
intrusive interference signal is present, it is analyzed 
to determine its frequency. An alternate spectral 
region is next identified that is free of interf erers . 
Lastly, communications are adjusted by re- tuning or 
15 modifying data rates to allow continued communications 
without intrusive interference. 

In another embodiment of the present invention, a 
modem capable of non-intrusively detecting an unwanted 
interferer while a receiving of a data signal is 
20 provided. The modem includes a transmitter for 
transmitting information to a second modem and a 
receiver for receiving the data signal from the second 
modem. A memory is connected to the receiver for 
storing statistical information associated with said 
25 data signal, including data sample values in each 

received symbol cluster in the data signal's amplitude 
phase constellation. A processor measures constellation 
statistics from the statistical information, determining 
a carrier-to-noise ratio and bit error rate of said data 
30 signal, determining whether the bit error rate exceeds 
an acceptable limit by testing whether symbol cluster 
changes from a Rayleigh distribution to a Ricean 
distribution, and when the bit error rate exceeds said 
acceptable limit, interupting the received data signal 
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to spectrum analyze a transmission band to find an 
interference free region of the transmission band and 
communicating the interference free region to the second 
modem for adjustment of transmission parameters. 

5 

Brief Description of Drawings 

FIGs. 1A and IB are an illustration of the effect of a 

sinusoidal tone interferer upon the symbol clusters 

in a quadrature amplitude modulation (QAM) 

10 constellation. 

FIG. 2 is a block diagram of a communication system 

using the present invention. 

FIG-. 3 is a flow chart of the processing step of 

measuring the constellation statistics. 

15 FIG. 4 is a flow chart of the processing step of 

determining the carrier-to-noise ratio (CNR) and 

bit error rate (BER) . 

FIG. 5 is a graph of the standard deviation of 

determined CNR versus size of the input data set 

20 for a 16-QAM signal with a CNR of 12.8 dB . 

FIG. 6 is a graph of the standard deviation of 

determined CNR versus size of the input data set 

for a 64-QAM signal with CNR of 19.1 dB . 

FIG. 7 is a graph of the ratio of the mean to standard 
25 deviation, pi(z), versus z for a Rice distribution. 

FIG. 8 is a graph of the skewness, p 3 (z), versus z for a 

Rice distribution. 
FIG. 9 is a graph of the kurtosis, p 4 (z), versus z for a 

Rice distribution. 
30 FIG. 10 is a graph depicting cluster skewness as a 

function of CNR and interferer amplitude versus 

inteferer power. 
FIG. 11 is a graph depicting cluster skewness as a 

function of CNR and interferer amplitude versus 
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normalized interferer amplitude. 

Best Mode for Carrying Out the Invention 

The present invention provides an improved method 
5 for recognizing the presence of a CW interferer in the 
demodulated output of a CAP or QAM digital data 
receiver. The method relies on the fact that the vector 
combination of the interferer with a CAP/QAM symbol 
cluster changes the radial distribution of sample 
10 locations in the cluster from a Rayleigh distribution in 
the presence of only Gaussian noise to a Ricean 
distribution in the presence of such an interferer. The 
invention uses statistical analysis of the positions of 
a set of samples in the demodulated symbol constellation 
15 to decide whether the distribution is sufficiently 

shifted towards Ricean to declare that an interferer is 
present . 

Three different measures of the shape of a 
probability distribution function were tested: the ratio 
20 of the mean to standard deviation, the skewness or third 
central moment, and the kurtosis or fourth central 
moment. Of these measures, the skewness was found to be 
best correlated with the received carrier-to- 
interference ratio (CIR) , and the kurtosis was also 
25 found to provide a useful degree of correlation. By 
calculating either the average skewness or average 
kurtosis over all clusters in the constellation, and 
testing this average skewness or average kurtosis 
against a threshold, the present invention determines 
30 that an interferer is present if the average skewness or 
average kurtosis is below the threshold, and that no 
interferer is present if the average skewness or average 
kurtosis is above the threshold. 
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While a fixed threshold is possible, it has been 
found that at high CNRs the interference levels 
(measured as a fraction of the noise margin) that can be 
positively recognized in this way are higher than the 
5 interference levels that can be recognized at low CNRs. 
In high noise (low CNR) situations there is not much 
difference in the skewness or kurtosis of an 
interference-free signal and one that has a moderate 
level interferer that could otherwise be detected in a 

10 high CNR environment. With a small set of sample data, 
the statistical fluctuations in the computed skewness or 
kurtosis can be larger than the significant differences 
in skewness. However, with a large enough sample set, 
even small differences in skewness or kurtosis and con- 

15 sequent CIR can be detected, even at low CNR, providing 
that the threshold criterion is made a function of the 
CNR (which is routinely measured) . Generally, skewness 
has been found to be more stable with respect to 
sampling fluctuations than is kurtosis. The main 

20 limitations will be the speed with which the data 
samples can be transferred to the data processing 
computer and the processing speed of the computer. For 
a cost-limited unit with limited computing resources, 
such as a VDSL subscriber modem, there may be design 

25 limitations on the number of data samples or the 
sampling rate that can be handled. 

According to the present invention, the receiver 
accumulates certain statistics regarding the data sample 
values in each received symbol cluster in the signal's 

30 amplitude-phase constellation. These statistics are 
used to calculate the CNR (actually, the carrier to 
noise-plus-interference ratio) and the BER for the 
received signal constellation. If the BER falls below 
an acceptable limit, then the average cluster skewness 
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and/or kurtosis is also computed for the constellation. 
A table lookup provides the threshold skewness and 
kurtosis values at the calculated CNR below which a 
calculated skewness or kurtosis would indicate the 
5 presence of an interferer of amplitude sufficient to 
significantly affect the measured CNR. The calculated 
skewness and/or kurtosis are then compared to the 
selected threshold skewness and/or kurtosis, and if the 
calculated skewness and/or kurtosis are below the 
10 threshold, then the presence of an interferer is 
declared. 

The frequency of the interferer is then determined 
by spectral analysis of the received transmission band 
for the signal of interest using, most commonly, an FFT. 

15 If the interferer is strong (small CIR) , the spectral 

line of the interferer may show up above the spectrum of 
the SOI. However, if the amplitude of the interferer is 
approximately the same as that of the SOI, it will 
generally be necessary to interrupt transmission of the 

20 SOI long enough to collect interference data to spectrum 
analyze. In some cases, it is possible to infer the 
frequency of the interferer if there can only be one 
such frequency. For example, in VDSL transmission there 
may be only one allocated amateur radio band within the 

25 spectral band used for a particular VDSL transmission in 
which interference is detected. This will then be the 
interferer frequency. It may then be possible and 
acceptable to avoid the interferer by reducing the 
transmission bandwidth (and the corresponding data rate) 

30 to the segment of the original transmission band on one 
side of the interferer, so that the interferer is then 
outside the new VDSL transmission passband. 
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If reducing the data rate is not an acceptable 
alternative, then it will be necessary to find another 
passband within the available VDSL spectral range that 
is free of interferers. This requires spectrum 
5 analyzing the entire VDSL spectral range, not just that 
passband in current use. This must be done before any 
bandpass filtering in the receiver, and will again 
require interrupting the VDSL data transmission long 
enough to collect data for spectral analysis of the 

10 interference environment. The advantage of the present 
invention is that this interruption of the SOI data 
transmission need only be done when it has been 
determined from constellation analysis that there is an 
interferer present which is causing unacceptable 

15 degradation of the BER. If the BER is so degraded, then 
it does not matter that data transmission is 
interrupted, since the received data already has too 
many errors. It is better to interrupt transmission and 
find an acceptable transmission band so that suitably 

20 error-free transmission can be restored as quickly as 
possible. 

The present invention detects when an interferer is 
present that is degrading the BER of the SOI. The 
invention then provides means for determining a 

25 transmission band that will be free of degrading 

interferers and then retuning the transmitter to the new 
transmission band. 

Referring now to the figures, and for the present 
to Fig. 2, therein is shown a block diagram of the 

30 present invention. In the transmitter 16, the input bit 
stream 2 0 enters the symbol encoder 22 where the input 
bit stream 2 0 is encoded into QAM /CAP modulation. The 
resulting symbol sequence has in bits per symbol in a 
constellation of M symbols. Generally, m=log 2 M. In the 
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modulator 24, the symbol sequence out of the symbol 
encoder 22 is modulated onto in-phase and quadrature 
rails at a carrier- frequency derived from carrier 
oscillator 26. The modulated signal is then boosted in 
5 power by the power amplifier 28 and the signal of 
interest (SOI) is transmitted over the forward 
transmission channel 30 to the receiver 18. Along the 
forward transmission channel 30, external radio- 
frequency interference 32 is added to the transmitted 
10 SOI, the addition being modeled by the summation 34. 

The combined signal of SOI plus RFI enters the 
receiver 18 where it is first demodulated in the 
demodulator 36. This includes locking to the carrier 
frequency and symbol clock rate, separating the in-phase 
15 and quadrature components, downconverting from the 

carrier frequency to baseband, and sampling the signal 
at the symbol clock rate. The output of the demodulator 
3 6 is a sequence of in-phase and quadrature component 
values for the sum of SOI plus noise plus interference 
20 at each symbol sample. The output of the demodulator 3 6 
is passed through the equalizer 38, which compensates 
for and removes the effects of intersymbol interference 
and changes in sample values caused by channel 
distortion, filtering, and multipath propagation (e.g., 
25 extraneous paths involving bridged taps, i.e., open 
circuited twisted pairs connected in shunt with the 
main-path twisted pair on a VDSL line) . The equalized 
signal goes to slicer 40 which tests the in-phase and 
quadrature values against decision boundaries and makes 
30 a decision as to the symbol state for each equalized 
sample. The symbol stream 42 out of the slicer then 
goes to the symbol decoder 44. The symbol decoder 44 
converts each symbol to the binary bit sequence 
corresponding to the input symbol state and outputs a 
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binary bit stream 46 representing the information 
conveyed by the system. The bit stream 46 goes to an 
interpreter (not shown) which separates the payload 
information from the overhead control bits. 

5 The signal at the output of the equalizer 3 8 is 

sampled by the left hand pole 66 of double-pole switch 
64, and the signal at the output of the slicer 40 is 
sampled by_ the right hand pole 68 of double-pole switch 
64. Both poles of the switch 64 are sampled in 

10 synchronism under control of receiver controller 62 . 
Both sets of sample data are stored in data storage 
buffer 72 within computer 76, with the output samples of 
equalizer 3 8 being delayed so that they are aligned with 
the corresponding output samples of the slicer 4 0 for 

15 the same symbols in parallel columns of the buffer 72 . 

These sample times are determined by control bits in the 
output bit stream 46 and by constraints imposed by the 
data collection capabilities and processing rates of the 
signal processing path within computer 76. For example, 

20 the switch 64 might be commanded to sample only every n- 
th symbol out of the equalizer 3 8 and its corresponding 
sliced value out of slicer 40. In addition, as will be 
described subsequently, switch 64 will normally be 
commanded to collect blocks of sample data only at 

25 specific times as will be described. 

The data storage buffer 72 consists of two blocks 
of random access memory (RAM) within computer 76. The 
two blocks of RAM are "pingponged" , that is, while data 
is being loaded into one block, the data in the other 

30 block is being processed. When the data in the one 

block has been processed, the contents of that block are 
overwritten with new data and the data in the other 
block are then processed. As is known by those skilled 
in the art, other types of storage can be used in place 
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of the RAM, including, for example, EPROM, SDRAM, hard 
disk, etc. 

At system startup, data are initially transferred 
over the forward transmission channel 3 0 using a 
5 relatively low order constellation, say 16-QAM, with a 
relatively low symbol rate and with a bandwidth and 
carrier frequency selected to have a high probability of 
being free of RFI. At startup, the receiver 18 is 
automatically set by the receiver controller 62 to tune 
10 to the standard startup carrier frequency, symbol rate, 
and constellation size. Thereafter, prior to any time 
that the transmitter 16 plans to change either the 
carrier frequency, symbol rate or constellation size, 
the transmitter 16 sends a transmission change 
15 notification to the receiver 18. The receiver controller 
62 extracts the new transmission parameters from the 
output bit stream 4 6 and sends an acknowledgment over 
reverse control channel 58 to the transmitter controller 
60. Only when the transmitter controller 60 has 
20 received this acknowledgment does it command the 

transmitter 16 to change the transmission parameters. 
At this point the receiver controller 62 has commanded 
the receiver 18 and the computer 7 6 to expect the new 
transmission parameters. Knowledge of the new carrier 
25 frequency and symbol rate permits the demodulator 3 6 to 
rapidly acquire lock to these new parameters, while the 
slicer 40 decision boundaries are reset to correspond to 
the new constellation. 

Once initial receiver synchronization has been 
30 received and acknowledged to the transmitter controller 
60 over the reverse control channel 58, the transmitter 
controller 60 informs the receiver controller 62 over 
the forward transmission channel 3 0 and then turns off 
power amplifier 28 to interrupt the transmitted signal: 
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The receiver carrier frequency and symbol rate phase 
locked loops in demodulator 3 6 continue to flywheel at 
the previously acquired carrier frequency and symbol 
rate during this interruption. Also during the 
5 interruption, switch 70 is closed and switch 64 is 
opened, so that the radio frequency interference 
environment is sampled and digitized by analog-to- 
digital converter 78. Computer 76 then collects 
digitized time samples of the interference environment 
10 which it spectrum analyzes using an FFT 74. The 

resulting interference environment spectrum is trans- 
mitted over reverse control channel 58 to transmitter 
controller 60, which selects a new carrier frequency, 
symbol rate and constellation size that are adequate to 
15 carry the desired payload data rate while avoiding any 
RFI that may be present in the environment. After 
notifying the receiver controller 62 of the new 
parameters and receiving an acknowledgment, the 
transmitter controller 60 resets the parameters in the 
20 symbol encoder 22 and carrier oscillator 26, and turns 
back on the power amplifier 28 to begin data 
transmission. 

The first processing step on any block of input 
data out of data storage buffer 72 is to measure the 
25 constellation statistics 48. The flow chart for this 
processing step is shown in Figure 3. The computer 76 
receives information in step 3 02 from the receiver 
controller 62 as to the size of the symbol 
constellation. A block of data storage is then 
30 allocated to each symbol cell at step 304, and the 

symbol data out of the equalizer 3 8 are sorted into the 
appropriate symbol cell locations at step 306 according 
to the corresponding cell identifiers that come from the 
output of the slicer 40. After the presliced symbol 
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data is sorted into bins by cell, the number of symbols 
in each cell are counted (step 308), and the sample mean 
values and sample variances of the data in the in-phase 
(x) and quadrature (y) directions are calculated for 

5 each bin at steps 310 and 312. This may be done for 
each symbol cell in the constellation or, for large 
constellations such as 256-QAM, for just the cells in 
one quadrant. With M cells and 1^ sample symbols in 
each cell, the sample means and sample variances are 

10 given by: 

(1) Sample mean x-position of sample set in cell m: 

1 "™ 
IN mn=i 

(2) Sample mean y-position of sample set in cell m: 

,N m n =l 

15 (3) Sample variance on x-positions of sample set in 

cell m: 

r& , _ i 



—L—'v 2 _- 

N m -iLs Xmn " Xm J' 



(4) Sample variance on y-positions of sample set in 
cell m: 

1 THn, 1 

20 £y.m = "j^— j ^2,YFnn- YmJ • 

The numbers of samples N m in each cell are counted 
at step 3082, since these numbers are needed in the 
above calculations . 

The next processing step (see block 50 of FIG. 2) 
25 is the calculation of the received signal carrier-to- 
noise ratio (CNR) and probability of error, or bit error 
rate (BER) . The CNR is determined using the method of 
symbol cluster size analysis, wherein the CNR is 
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determined as a function of the constellation closure 
ajd, where a n is the RMS noise voltage (equal to the RMS 
symbol cluster size) and d is the noise margin, i.e., 
the minimum (perpendicular) distance, in volts, from the 

5 ideal symbol point to the nearest decision boundary. The 
quantity is called the cluster variance. This method 
is usable for high enough CNRs that each symbol cluster 
is entirely within its decision cell. A flow chart for 
the procedure is shown in Figure 4. First at step 402, 

10 the RMS cluster size a n is calculated using the 
following equations : 

M 

N=lN m , 
rn=i 

M 



15 



varx= 7^£ (Nm - 1)s; 



1 M 

vary -TTT^£, (H "- 1)s 5. m - 



jvarX + 
r n = y J 



vary 



Then at step 404 the reciprocal of the constellation 
closure d/a n is calculated. The CNR is then obtained at 
step 406 as 



CNR(dB) =201og 10 

20 where L is given for each constellation modulation in 
the following table: 
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Modulation 


L 


a 


b 


QPSK 


2 


2 


1/2 


16 -QAM 


10 


3 


3/4 


32-Cross QAM 


20 


13/4 


11/13 


64 -QAM I 


42 


7/2 


7/8 


128-Cross QAM 


82 


29/8 


53/58 


256 QAM 


170 


15/4 


15/16 


9-QPR 


5.3333 


8/3 


2/3 


25-QPR 


16 


16/5 


4/5 


49-QPR 


32 


24/7 


6/7 


81-QPR 


53 .3333 


32/9 


8/9 



TABLE 1 



The symbol error rate (SER) is obtained at step 408 
5 by first calculating the quantity Q: 

where erfc denotes the complementary error function. A 
lookup table for ef rc ( ) can be stored in the memory of 
computer 7 6 and interpolated for any value of the 

10 argument. The SER is then calculated using the following 
equation with the appropriate coefficients, as listed in 
the above table: 

SER = agi-bQ. 

As an alternative to using a lookup table for 

15 erfc, the function Q(u) can be evaluated by the 
polynomial approximation 

Qto> = ^= e u " /2 (qt + qt 2 + qt 3 + Qt 4 + qt 5 ), 

where 

1 

t = - , 

1+pu 
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\i=cfa n , and p and C\, . . . , c 5 are constants given by 

p = 0.2316419 

ci = 0.319381530 

c 2 = -0.356563782 
5 c 3 = 1.781477937 

c 4 = -1.821255978 

c 5 = 1.330274429 
The inherent accuracy of this approximation is better 
than 1% down to an error probability of 10~ 15 , and 10% 
10 down to 
lO" 148 . 

Next, the BER is calculated at step 410. For full 
response QAM and cross-QAM constellations having M 
symbol cells, the BER is related to the SER by 
1 5 BER = SER / bits per symbol , 

where 

bits per symbol = log 2 M. 

To use the constellation closure to determine CNR, 
and from it the SER and BER, it is necessary first to be 
20 able to resolve the individual symbol clusters, and 
second, to be able to determine the cluster variance 
with reasonable accuracy. It is not possible to resolve 
the individual clusters when the RMS cluster size cr n is 

close to the noise margin d, that is when ajd^l. In 
25 practice, it is not possible to resolve the individual 
symbol clusters when a T Jd>K, where K is some number in 
the range 0.5 < K < 1.0. When cr n /d=K, then 

(CNR) dB = 101og l0 (iy2K 2 ). 

This gives the limiting minimum value of CNR below which 
30 it is not possible to meaningfully determine the CNR for 
L corresponding to a particular constellation size and 
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for an assumed value of K. These limiting values of CNR 
are listed in the following table 2: 



\A r\H i i lit nn 




CNR(dB) 


K = 0.25 


K = 0.5 


K = 0.75 


K= 1.0 


QrbK. 


2 




o.u 


Z.j 


U.O 


16-QAM 


10 


19.0 


13.0 


9.5 


7.0 


32-Cross QAM 


20 


22.0 


16.0 


12.5 


10.0 


64-QAM 


42 


25.2 


19.2 


15.7 


13.2 


128-Cross QAM 


82 


28.2 


22.1 


18.6 


16.1 


256 QAM 


170 


31.3 


25.3 


21.8 


19.3 


9-QPR 


5.3333 


16.3 


10.7 


6.8 


4.3 


25-QPR 


16 


21.0 


15.0 


11.5 


9.0 


49-QPR 


32 


24.0 


18.0 


14.5 


12.0 


81-QPR 


53.3333 


26.3 


20.3 


16.8 


14.3 



5 TABLE 2 

To achieve reasonable accuracy in determining the 
CNR from cluster variance analysis, a large enough 
sample symbol set must be used to reduce the expected 

10 uncertainty, or standard deviation, of the calculated 
CNR to a reasonably small value. Specifying or 
measuring the CNR to an accuracy of much less than 0.1 
dB is generally not meaningful. Therefore, a large 
enough sample set must be used to assure that the 

15 uncertainty in the calculated CNR is no larger than 0.1 
dB. Some Monte Carlo simulations have been done with 
data sets of various sizes to estimate the standard 
deviation of the calculated CNR for various size 
constellations for a CNR that is theoretically adequate 



WO 00/07302 



19 



PCT/US99/16864 



to achieve a BER of lxlCT 7 (the maximum allowed error, 
rate for a VSDL modem) . 

FIG. 5 shows the standard deviation of calculated 
CNR versus size of the input data set for a 16-QAM 
5 signal with a CNR of 12,8 dB (needed to achieve a BER of 

lxlO" 7 ), and Figure 6 shows the standard deviation of 
calculated CNR versus size of the input data set for a 
64-QAM signal with a CNR of 19.1 dB (needed to achieve a 

BER of lxlO" 7 ). In each case, the standard deviation was 
10 calculated using a set of ten simulations at each input 
data set size, which accounts for the still jagged 
nature of each graph. Figure 5 shows that for 16-QAM at 
least 1500 symbol samples are needed, and for 64-QAM at 
least 3 000 symbol samples are needed to achieve a 
15 standard deviation of 0 . 1 dB in CNR at the critical CNR 
level. A similar analysis was done for 256-QAM, which 

needs a CNR of 25.15 dB to achieve a BER of lxlO" 7 ; it 
was found that 50,000 samples were needed to achieve a 
standard deviation of 0 . 4 dB in the calculated CNR, and 
20 it is estimated that 0.2-0.3 million samples would be 

needed to achieve a BER of lxlO -7 . Consequently, there 
will be a longer time needed between updates of the 
estimates of CNR and BER for 2 56-QAM than for smaller 
constellations. 

25 For the step of testing the BER (52 in FIG. 2), the 

BER that is calculated using the above described 
procedure is compared to the maximum allowed threshold 
BER. If the calculated BER is at or below the threshold, 
no change is made to the signal transmission parameters, 

30 and a new set of data is accumulated to update the 

estimated BER and repeat the test. If the calculated 
BER is above threshold, this suggests that there may be 
an interference source that is raising the error rate. 
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To check this, the computer 7 6 goes to the step of 
calculating the average cluster skewness and/or kurtosis 
54, rather than updating the estimated BER. It must be 
noted that the sample data in either RAM block of the 
5 data storage buffer 72 must not be overwritten with new 
data until the testing of the BER 52 is completed, since 
the data stored in that RAM block will be needed in the 
subsequent calculation of skewness and/or kurtosis 54 if 
the BER test 52 indicates that the BER has risen to an. 

10 unacceptably high value. To assure this, whenever the 
BER test 52 registers above threshold, the computer 76 
sends a signal informing of this condition to the 
receiver controller 62. The receiver controller 62 then 
opens both poles of switch 64 to prevent new data from 

15 entering into data storage buffer 72. 

When there is no interference present, the x and y 
components of samples in cluster m are each Gaussian 
distributed due to noise. If there are no distortions, 
x and y correspond to the I and Q components . 

20 Furthermore, if there is no gain mismatch, the x and y 
distributions will have the same variance. In that 
case, the radial distance to the samples from the center 
of the cluster will have a Rayleigh distribution. 

When there is an interferer present (assumed to be 

25 a sinusoidal tone at other than the digital carrier 
frequency) , the sinusoidal interferer causes each 
cluster to be ring shaped due to the addition of a 
randomly rotated interference vector of constant 
magnitude to each sample, as shown in Figure 1. If 

30 there were no noise on the system, the interferer would 
transform each ideal constellation point to a sharp 
circle of zero thickness. On the other hand, if the 
interferer had finite amplitude modulation, the locus of 
the interference vector would have a finite thickness. 
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This thickness would add to the width of the noisy, but 
interference-free cluster to result in an interference- 
plus -noisy- symbol cluster ring thicker than the ring 
that would result from a pure tonal interferer. 

5 When an interferer is present, the radial distance 

to the samples from the center of the cluster will have 
a Rice distribution. Thus, to automatically determine 
whether or not an interferer is present and affecting 
the digital modem performance, it is necessary to 

10 automatically distinguish between a Rayleigh and a Rice 
distribution, with the same peak radial distance on the 
Rice distributions for all symbol clusters. For the 
Rayleigh distribution, the probability density function 
is given by 

lWr)=-^e' J ' 2 *\ r>0, 
15 a 

where 

r = VOW-SrrOMymk- ?mf 
is the radial distance from the centroid fovYm) of 
cluster m to the k-th sample in the cluster, located at 
20 (^nkYmk)- Also the standard deviations of the cluster 
distribution in each coordinate are assumed equal, so 
that o- = cr v =a u . Using the same notation, the Rice 
distribution is given by 

IWtt-^e lo[p-J> '>0, 

25 where I 0 (x) is the zero-th order modified Bessel function 
of the first kind, and a the amplitude of the 
interfering signal. Here, a is the amplitude of the 
interferer. Note that when a = 0, I 0 (0) = 1, so that 
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To decide which type of distribution is present, 
and thus whether or not an interferer is present, 
various moments and functions of moments may be 
calculated for the sample distribution for each cluster, 
5 and these moment functions may be compared with the 
theoretical values for the Rayleigh and Rice 
distributions. This will work as long as the 
constellation distortion is not severe, so that the 
Rayleigh and Rice distributions are approximately valid 
10 for their appropriate situations. 

The n-th moment of a probability distribution p x (x), 
in the random variable X, is here denoted 

^;=E[x n ], (4) 

n a = E\(x- M ;f]. (5) 

15 /i 2 is the variance of the distribution. One measure that 
is useful for characterizing probability distributions 
is the ratio of the mean to the standard deviation, 

Px=ti/4ih. (6) 
Two other useful measures are the skewness, 

20 p y =nJ(M 2 )\ (7) 

and the kurtosis, 

P.=vJv; ( 8) 

(see N. L . Johnson and S. Kotz, Distributions in 
25 Statistics, Vol. 1, Discrete Distributions , Houghton 
Mifflin Company, 1969, p. 18) . The skewness is a 
measure of the asymmetry of the distribution. Positive 
values of p 3 indicate that the distribution is "skewed to 
the right" so that the right tail is in a certain sense 
30 heavier than the left. Similarly, negative values 

indicate that the left tail is heavier than the right. 
Kurtosis is a measure of the peakedness of the 
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distribution. For a normal or Gaussian distribution, p 4 
= 3. For p 4 < 3, the distribution is less peaked than a 

Gaussian; for p 4 > 3, it is more peaked. 

Equations (6) to (8) are expressed in terms of 
5 central moments, except for ji{ since there is no fi[. Note 
that the shape factors p x§ p 3 , and p 4 can all be expressed 

in the form p n = /i n (ju 2 ) 2 . However, p 2 is meaningless since 

Pi - 1 • 

The central moments can be expressed in terms of 
10 the moments (4) by 

^ =/l / 3 -3Ai 2 X+2/i I /3 (9) 



For a Rayleigh distribution, as in (1) , the 
moments (4) are given by 

^ = (2<r j rr(i+|]. (io) 

where T(k) is the Gamma function, with properties for 
k > 0, 

JkT(k) in general 
integers 



15 



20 



r( k + i)=j^ (lla) 



Inserting Equations (11) into (10) gives 

K = jfor, (12a) 
H',=2<?, (12b) 
^ = 3^<r\ (12c) 

^;=8ct 4 . <i2d) 

25 Then, inserting Equations (9) and (12) into Equations 
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(6) to (8) and simplifying gives the shape factors for 
the Rayleigh distribution as 



P '~ = -/^~ = 1 - 9131 ' (13a) 

P,^ = ^# = 0.6311, (13b, 
32- 3k 2 

p 4R , <lMn =— r =3.2451. (13c) 

5 R *' toRh (4-;r) 

For the Rice distribution in Equation (3), which is 
equivalent to a noncentral chi-square distribution with 
two degrees of freedom, the moments (4) are given by 

10 where ,F;(a,/?;X) is the confluent hypergeometric function, 
defined as 

,F(«^x)_X r(a)r ^ + k)k! 0*0,-1,-2,.. . 

Comparison of Equations (10) and (14) shows that the 
moments of the Rice distribution are related to those of 
15 the Rayleigh distribution by 



a' = u' e <f ' 2a F 



— + l;l; — , 
2 2cr J 



(15) 



Furthermore, defining 



CH 

Z=— 7, (16) 

2cr~ 



we have that as z-> 0, e~ z ^ 1 and ,F;(a; l;z) -> 1 for arbitrary 
20 a, so that and Equation (14) reduces to 

Equation (10) . Note that z is the interference-to-noise 
ratio (INRJ . 

Using published transformations of the confluent 
hypergeometric function, the moments of the Rice 
25 distribution may be expressed as 
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*<* = "/f e_Z L Fl(-°-5:l;z) + 2z ,F, (0.5;1; 2 )> = Q(z)<7, ( 17a > 
^R, ce = 2 < 1 + z)e z ,F,(l;l;z)cr 2 = Q(z)<r 2 , (17b) 



^ Wce = ^fe z [2(2 + Z) , F,(-0.5;l;z) + (4z 2 + 8z- 1), F, (0.5; l;Z)]r 3 =C3(Z)a ? , 
5 (17c) 

^4 Wce =8© 2 | F,(3;l;z)<r 4 = 4(z 2 + 4Z+ 2)x 4 = C 4 (Z)<x 4 , (!7d) 
where ,F,(-0.5;1;Z) , ,F,(0.5;1;Z) and ,F,(1;1;Z) are all tabulated 
functions (see M. Abramowitz and I. Stegun, Handbook of 
Mathematical Functions, Dover Publications, Inc., 1965, 
10 pp. 504-534) . Inserting these expressions into 

Equations (9) gives the central moments for the Rice 
distribution as 

H 2 (z)=\Q(Z)-[q{Z)f}o\ (18a) 
H 3 (z) = -fc(z)-3Q(Z)q(z) + 2[q(Z)f]o-\ (18b) 

l5 /i 4 (z) = -jc 4 (z)-4q(z)q(z)+6Q(z)[q(z)] 2 -3[q(z)] 4 ]a 4 . (isc) 

Then, inserting Equations (17a) and (18) into Equations 
(6) to (8) gives the probability distribution shape 
factors p„. The results are 

p,(z)=-=HE_, (19a) 

VQ(z)-[q(z)]- 



20 



Q(z)-3Q(z)q ( z)+2q(z)- 
f;(z)-[q ( z)]-] 

(;:) c 4 (z)-4q(z)q ( z ) + 6C ?( z)[q ( z ) f -3[q ( z)] 4 
^ ( z)-[q(z)]^ 

Equations (19) are plotted in Figures 7 to 9 . 
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The Rice distribution is known to approach the 
Rayleigh distribution as z->0, and to approach a Gaussian 
distribution for large z. Note in Figures 7 to 9 that 
the curves for p n have the values given in Equations (13) 
5 for the Rayleigh distribution at z = 0 . For large values 
of z, the kurtosis, p 4 , approaches the limiting value of 
3, as for a Gaussian distribution. Similarly, the 
skewness, p 3 , approaches 0, indicating that the 
distribution is becoming more symmetric, like a Gaussian 

10 distribution. For small z, p 3 is large, indicating that 
the distribution is highly asymmetric. 

It appears from Figures 7 to 9 that a reasonable 
test as to whether a CW interferer is present in the 
data, and corrupting the constellation, is to calculate 

15 the ratios p, , P 3 , and p 4 , for the observed data averaged 
over all symbol clusters, and compare the values with 
Figures 7 to 9 . If p^HO, p 3 <0.50, and p 4 <3.10, one can 
be reasonably certain that the radial probability 
distribution of points in the clusters is not Rayleigh, 

20 and that there is probably an interferer present. Based 
upon the curves in these figures, if a single parameter 
were to be used to make this decision, it should be the 
skewness, p 3 , which exhibits a rapid and monotonic 
departure from the Rayleigh value at z = 0 . However, 

25 the kurtosis, p 4 , also exibits a useful distinction 
between large values for small z (small or no 
inter ferers) and smaller values for large z (large 
inter ferers) , although the kurtosis function is not 
monotonic . 

30 To check the validity of this approach, a number of 

simulations were run for a range of input CNRs and a 
range of interferer amplitudes, where the amplitude a 
was normalized to the noise margin d, and the normalized 
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interferer amplitude was varied over the range 0.000001 
< a/d < 0.5. Although most runs were done with an 
amplitude spread of 0 on the interferer (no amplitude 
modulation) , some runs were also done with a uniformly 
5 distributed amplitude spread of 0.1 to represent a 
modest amplitude distribution. Similar results were 
obtained with and without the amplitude modulation. The 
results of these simulations showed a strong correlation 
between the cluster skewness p 3 and the calculated 

10 carrier-to-interference ratio (CIR) , and a somewhat 

weaker correlation between the cluster kurtosis p 4 and 
the CIR. The distribution shape parameter p, showed a 
much poorer correlation with the CIR. This prediction 
that the cluster skewness is the best single predictor 

15 of the presence of an interferer was subsequently- 



had the capability to make the various measurements and 
tests. Consequently, although both the skewness and 
kurtosis are shown in Step 54 of the procedure in Figure 

20 

2, the skewness is the preferred parameter if just one 
is to be calculated. 

To calculate the skewness in Step 54, the symbol 
coordinates for each sample symbol in each cell m are 
adjusted by subtracting the sample mean coordinates 3^ 
25 and y m to form the coordinate differences from the 
cluster center, and y mk -? m . Then the radial 

distance of each sample from its cluster center, 



is calculated. If there is no distortion of the 
30 constellation, so that the expected values of the sample 
mean coordinates 3^ and y m coincide with the ideal 
locations of the cluster centers, then the sample means 



confirmed by 



tests on a developmental VDSL modem that 




(2 
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in Equation 20 may be replaced by the ideal cluster 
coordinates. This also eliminates statistical 
fluctuations in the calculated cluster mean coordinates. 
Next, within Step 54, the sample moments 

5 ^ = Tak\^^ t ^ n=l,...,4, 

are calculated, where N m is the number of data symbol 
samples in the bin for symbol cell m, accumulated 
previously in Step 48, and N is the total number of data 
symbol samples in all bins. Analogous to Equations (9), 
10 the second and third central moments are then calculated 
by 

S3 = s^-3s;s[ + 2s,' 3 
The cluster skewness is then calculated by 

15 ft = s$/(s 2 ^. 

Similarly, to calculate the kurtosis in Step 54, the 
fourth central moment is calculated by 

s 4 =S^-4^s ! , + 6s;s 1 /2 -3s[ 4 , 
and the cluster kurtosis is then calculated by 
20 p 4 = Sj/s 2 2 . 

To determine whether an interferer is present, the 
calculated skewness and/or kurtosis is compared, in 
Step 56, with a threshold level determined from Figure 8 

25 or 9, respectively. As seen from Figure 8, a skewness 
less than about 0.50 indicates that the cluster radial 
distribution is Ricean distributed, indicating that an 
interferer is present, while a skewness greater than 
this value indicates a Rayleigh distribution, indicating 

30 that no interferer is present. Similarly, from Figure 9 
it is seen that a kurtosis less than about 3.15 
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indicates that the cluster radial distribution is Ricean 
distributed, indicating that an interferer is present, 
while a kurtosis greater than this value indicates a 
Rayleigh distribution, indicating that no interferer is 
5 present . In experiments on the developmental VDSL 

modem, a skewness threshold of 0.3 to 0.5 and a kurtosis 
threshold of 3.12 to 3.15 have been found to be optimum. 
However, the lower the CNR of the input signal, i.e., 
the higher the noise level, the stronger the interferer 
10 must generally be to detect it above the noise. Of 
course, the smaller the interference-to-noise ratio 
(INR) is, the less proportionate effect the interferer 
has on the overall SER. When the cluster skewness is 
plotted as a function of both the CNR and the normalized 
15 interferer amplitude, as in Figure 10 (which was 

generated for a 4 -QAM signal), it is found that, except 
for some anomalous behavior in the vicinity of CNRs from 
15 to 20 dB (which is reproducible, but of undetermined 
cause) , a given skewness threshold level can only 
20 detect larger amplitude interferers as the CNR 

decreases. Similarly jagged and anomalous results have 
been found when using the cluster kurtosis as the 
interference detection measure, as shown in Figure 11. 
With either approach, neither the skewness nor the 
25 kurtosis can be used to determine the interferer 's 

amplitude. But either can be used to determine whether 
an interferer with amplitude significantly (more so for 
kurtosis than for skewness) higher than that 
corresponding to that of the test criterion is present 
30 (except at low input CNRs and high interferer amplitudes 
where the contours of skewness or kurtosis versus CNR 
are essentially vertical). Within these limits, for any 
CNR, a skewness and/or kurtosis threshold can be 
selected to determine whether an interferer with 
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amplitude greater than a specified level is present. 
For high noise and strong inter ferers, the combined 
effect is to expand the size of the symbol clusters so 
much that they overlap, at which point determination of 

5 CNR, SER, and cluster skewness and kurtosis from 

constellation analysis is no longer possible. Also, 
when there is much noise (low CNR) , the resulting 
statistical fluctuations cause some uncertainty in the 
calculated skewness and kurtosis, reducing the 

10 significance of the threshold decision. Nevertheless, 
at moderate and high CNRs, the method works quite well, 
and has been found valuable in detecting weak 
interf erers . 

Note from Figure 7 that the quantity Pi = sj/^s^ has a 

15 unique value for any value of the quantity z given in 
Equation (16), which is the INR . Therefore, if it has 
been concluded that an interf erer is present, the value 
of pi may be calculated, and the corresponding value of z 
may be looked up in a tabulated version of Figure 7 to 

20 obtain an estimate of the INR for diagnostic purposes. 

When the skewness and/or kurtosis test 56 has found 
the presence of an interf erer, the computer 7 6 sends a 
notification message via reverse control channel 58 to 
transmitter controller 60. Transmitter controller 60 

25 then interrupts the input bit stream 20 and turns off 
the power amplifier 28, so that the only input to the 
receiver 18 is the external interference 32. Meanwhile, 
the computer 76 also sends a notification message to 
receiver controller 62, which sets the carrier frequency 

30 arid symbol rate synchronizers to flywheel, opens both 
poles of switch 64, and closes switch 70, so that the 
received external interference 32 goes (through an 
amplifier, not shown) to the analog-to-digital 
conversion 78 and spectrum analysis 74. The receiver 
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controller 62 also sends an acknowledgment to the 
transmitter controller 60 that the receiver 18 is 
temporarily ceasing to process received digital data 
while the computer is performing the spectrum analysis. 

5 The transmitter controller 60 does not turn off the 
transmitter 16 until this acknowledgment is received. 
Simultaneously, the computer 76 switches its operation 
entirely over to spectrum analysis 74. The results of 
the spectrum analysis are sent via reverse control 

10 channel 58 to transmitter controller 60 which analyzes 
the received spectrum to find a carrier frequency at 
which the desired data rate may be transferred without 
encountering any interference. If no pass band of 
sufficient width is available, the transmitter 

15 controller 60 will switch transmission to a higher order 
constellation with a lower symbol rate, providing that 
such a combination is available that can provide the 
desired bit rate, and provided that the noise margin in 
the new constellation is adequate to provide the 

20 required BER. If these conditions cannot be met, the 
transmitter controller 60 will attempt to reduce the 
transmitted bit rate with the same or larger 
constellation in order to reduce the transmitted 
bandwidth, providing that the information to be 

25 transmitted will not be impaired by a reduced bit rate, 
as would some types of real-time transmissions. 

After selecting a new carrier frequency, 
constellation size, and/or symbol rate, the transmitter 
controller 60 inserts a control message into the input 

30 bit stream 20 to transmit over the forward transmission 
channel 30 to the receiver 18. The receiver controller 
62 extracts this message out of the output bit stream 46 
and then sets the receiver 18 to receive and demodulate 
a signal with these new parameters. The receiver 
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controller 62 also sends an acknowledgment of the 
receipt of this control message over the reverse control 
channel 58 to the transmitter controller 60. Only after 
receiving this acknowledgment does the transmitter 

5 controller 60 actually change the transmitter parameters 
to the new values . 

At this point, receiver controller 62 also opens 
switch 70, closes both poles of switch 64, and directs 
computer 76 to switch back to accumulating statistics on 

10 the new received symbol stream: calculating the CNR and 
BER; testing the BER; and, if necessary; calculating and 
testing the skewness and/or kurtosis on the new received 
constellation. 
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CLAIMS : 

1. In a communication apparatus, a method of detecting 
5 whether an interference signal is being received 

with a received signal of interest (SOI), said 
method comprising steps of: 

(a) accumulating statistical information 
associated with said SOI; 
10 (b) determining whether an intrusive 

interference signal is present with said 
SOI; 

(c) analyzing, when an intrusive interference 
signal is present, said intrusive 

15 interference signal to determine the 

frequency thereof ; 

(d) identifying a spectral region free of 
interferers; and 

(e) adjusting said communication apparatus to 
20 operate in said spectral region 

identified in step (d) . 

2. The method according to claim 1 wherein said steps 
(a) through (d) are non-intrusive to a continued 
25 reception of said SOI in said communication 

apparatus . 



3. The method according to claim 2 wherein said step 
30 (b) further comprises steps of: 

(f) determining a root mean square (RMS) 
cluster size of said intrusive 
interference s ignal ; 

(g) determining a reciprocal constellation 
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closure of said intrusive interference 
signal; and 
(h) determining a carrier-to-noise ratio 
(CNR) of said intrusive interference 
5 signal. 



10 



i. The method of claim 3 wherein said step (b) further 
comprises determining a symbol error rate of said 
intrusive interference signal. 

5. The method of claim 3 wherein said step (b) further 
comprises determining a bit error rate of said 
intrusive interference signal. 

15 6. The method of claim 3 wherein said step (b) further 
comprises determining an average cluster skewness ■ 
associated with said intrusive interference signal. 

7. The method of claim 3 wherein said step (b) further 
20 comprises determining an average cluster kurtosis 

associated with said intrusive interference signal. 



8. A transceiver for transmitting and receiving data 
signals, said transceiver being capable of non- 
25 intrusively detecting an interference signal to said 
data signals, said transceiver comprising: 

storage means for accumulating statistical 
information associated with said data signals; 

interference determining means for determining 
30 whether an unacceptable interference signal is present 
with said data signals; 

analysis means for spectrum analyzing said 
unacceptable interference signal and determining a 
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frequency thereof and further locating a spectral region 
free of interferers; and 

conveyance means for communicating to a 
transmitter said spectral region free of interferers for 
5 changing operation of the transceiver and transmitter to 
said spectral region free of interferers. 

9 . The transceiver of claim 8 wherein said 
interference signal is a continuous wave signal. 

10- The transceiver of claim 9 wherein said transceiver 
uses quadrature amplitude modulation. 

11. The transcevier of claim 9 wherein said transceiver 
15 uses carrierless AM/PM data communication methods. 

12. The transceiver according to claim 9 wherein said 
storage means accumulates both presliced quadrature 
amplitude coordinates and sliced symbol cell 

20 identification data. 

13 . The transceiver according to claim 12 wherein said 
analysis means processes a sequence of demodulated 
symbols to calculate the statistics of in-phase and 

25 quadrature mean coordinates and in-phase and quadrature 
sample variances for a subset of symbols in each 
decision cell from said sequence of demodulated symbols, 
said sliced data being used to assign said presliced 
data to a proper symbol cell. 

30 

14. The transceiver according to claim 13 wherein said 
analysis means further determines from said statistics a 
root mean square (RMS) cluster size for all the symbol 
clusters in a constellation, a constellation closure, a 
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resulting carrier-to-noise-ratio, a resulting symbol 
error rate, and a resulting bit error rate (BER) . 

15. The transceiver according to Claim 14 wherein said 
5 analysis means further tests the magnitude of said bit 

error rate to determine whether the quality of said 
demodulated signal is above an acceptable threshold and 
if not, then directing the commencement of further 
processing to determine whether the degradation of the 
10 quality of said demodulated signal is due to the 

presence of said continuous wave interfering signal. 

16. The Transceiver according to claim 14 wherein said 
analysis means determines an average cluster skewness 

15 for the symbol clusters in said constellation and tests 
the magnitude of said skewness to decide whether an 
interferer is most likely present, which would be the 
cause of said degradation of said BER. 

20 17 . The transceiver according to claim 14 wherein said 
analysis means determines an average cluster kurtosis 
for the symbol clusters in said constellation and 
testing the magnitude of said kurtosis to decide whether 
an interferer is most likely present, which would be the 

25 cause of said degradation of said BER. 

18. The transceiver according to claim 14, wherein 
said conveyance means transmits the results of said test 
of skewness back to said transmitter over a reverse 
30 channel, said transmitter temporarily turning off 
transmission while said receiver receives only the 
external interference signals in the RF environment, and 
spectrum analyzing said environment to determine the 
frequencies of said interference signals. 



WO 00/07302 



37 



PCT/US99/16864 



19. The transceiver according to claim 14 wherein said 
conveyance means transmits the results of said test of 
kurtosis back to said transmitter over a reverse 

5 channel, said transmitter temporarily turning off 
transmission while said receiver receives only the 
external interference signals in the RF environment, and 
spectrum analyzing said environment to determine the 
frequencies of said interference signals. 

10 

20. The transceiver according to claim 14 wherein said 
conveyance means transmits the results of said spectrum 
analysis back to said transmitter over a reverse 
channel, said transmitter determining new transmission 

15 parameters to avoid said interference signals, and said 
transmitter transmitting said new transmission 
parameters to said transceiver. 

21. The transceiver according to claim 20 wherein said 
20 transciever and said transmitter reconfigure to transmit 

and receive the desired data signals with said new 
transmission parameters and then resuming data 
transmission using said new transmission parameters. 

25 

22. The transceiver according to claim 20 wherein said 
new transmission parameters include a new carrier 
frequency, QAM modulation order, and data rate. 

30 23. The transceiver according to claim 9 wherein said, 
transceiver uses quadrature partial response modulation. 

24. A modem capable of non-intrusively detecting an 
unwanted interferer during a receiving of a data signal, 
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comprising : 

a transmitter for transmitting information to a 
communications apparatus; 

a receiver for receiving said data signal from said 
5 communication apparatus; 

a memory coupled to said receiver for storing 
statistical information associated with said data 
signal; and 

a processor for measuring constellation statistics 
10 from said statistical information, determining a 

carrier-to-noise ratio and bit error rate of said data 
signal, determining whether a continuous wave or 
narrowband interferer is present, determining whether 
said bit error rate exceeds an acceptable limit by 
15 testing whether a symbol cluster changes from a Rayleigh 
distribution to a Ricean distribution, and when said bit 
error rate exceeds said acceptable limit, interrupting 
said received data signal to spectrum analyze a 
transmission band for an interference free region of 
20 said transmission band and communicating said 

interference free region to said communications 
apparatus for adjustment. 

25 25. The modem according to claim 24 wherein said 

testing comprises calculating a skewness associated with 
said symbol cluster. 

26. The modem according to claim 24 wherein said 

30 testing comprises calculating a kurtosis associated with 
said symbol cluster. 

27. The modem according to claim 23 further comprising 
an analog-to-digital converter coupled for converting 
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said data signal into a digital signal. 

28. The modem according to claim 25 further comprising 
a spectrum f ast -Fourier-transf orm-analysis means for 
5 spectrum analyzing said analog signal. 
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